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On Collision Invariants for Linear Scattering 


By Laure Saint-Raymond and Mark Wilkinson 


Abstract 


In this article, we extend the result of Boltzmann [5] on characterisation of collision invariants from 
the case of hard disks to a class of two-dimensional compact, strictly-convex particles. 


1. Introduction 

Understanding the statistical behaviour of dynamical systems comprised of identical inter¬ 
acting particles has been a well-studied problem since the work of Boltzmann [5] in the kinetic 
theory of rarified gases. By studying the precise way in which particles scatter from each other 
following a collision, one is able to derive information about macroscopic properties of the system, 
such as the evolution of the local density of the gas or local propagation of heat. A great portion 
of the kinetic theory literature is devoted to the study of systems in which the identical particles 
are perfect spheres. However, it is a very natural question to understand in what ways the statisti¬ 
cal properties of systems of non-spherical particles differ from those composed of their perfectly 
spherical counterparts. 

In this article, we offer a preliminary contribution to the extension of the theory of the Boltz¬ 
mann equation from hard spheres to general hard particles. In the first part of this work, we study 
the physical dynamics of compact, strictly-convex bodies which do not interpenetrate. Moreover, 
we restrict our attention to systems of two identical particles, thereby considering binary particle 
interactions alone. The first important step in studying such systems is to construct suitable physi¬ 
cal boundary conditions for a dynamics (by means of scattering maps) when the two hard particles 
collide, in order that trajectories in phase space may be defined globally in fime. By ‘physical’ 
boundary conditions, we mean fhaf (i) fhe particles should not interpenetrate following collision, 
and (ii) there should also be conservation of total linear momentum, angular momentum and ki¬ 
netic energy of the two particles through any collision event. However, it is important to note here 
that, according to Wilkinson [17], it is not possible to construct a family of scattering matrices 
corresponding to the collision of two non-spherical particles which conserves their total linear mo¬ 
mentum, angular momentum and kinetic energy. Nevertheless, with the extension of Boltzmann’s 
equation to systems of non-spherical particles in mind, we construct and study families of scatter¬ 
ing matrices for two particle systems which conserve total linear momentum and kinetic energy of 
the colliding particles. 

The second and principal part of this paper is devoted to the important topic of collision 
invariants for non-spherical particle scattering in kinetic theory. To illustrate the importance of 
collision invariants, we turn very briefly fo fhe case of hard particles with spherical symmetry in 
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(which are hard disks in the case d = 2, but hard spheres in the case d = 3) and the classical 
Boltzmann equation. 


1.1. The Boltzmann Equation and Collision Invariants. It is well known that the Boltz¬ 
mann equation for the 1-particle density function / = f{x, v, t) given by 

^ + (V • V,)/ = C(/, /) for {X, V) e X (1) 

is a candidate PDE to describe the statistical properties of systems of N hard particles with spher¬ 
ical symmetry in the Boltzmann-Grad limit as ^ oo and e ^ 0 with Ne^~^ = 1, where £ > 0 
denotes the radius of any given particle. The unique family of scattering matrices {cr„)„ggrf-i which 
resolves a collision between two spherical particles, in such a way that properties (i) and (ii) above 
are satisfied, are the reflection matrices 

(Tn - 27 „ g) y„ € 0(2d), 

with where n € denotes the direction connecting the centres of mass of the 

two spheres at collision. The collision operator C{f,f) that appears in (1) is given by 

Cif, f) ^ f f |(v - v) • nl(f(x, v'„, t)f{x, %,t)- fix, V, t)fix, v, t)) dndv, (2) 
^ Jr”' Jsf' 

where the ‘post-collisionaT velocities [v^, v^] cr„[v, v] e R^^ are 

= V - [(v - v) • n]n and = v -i- [(v - v) • n]n. 


In order to derive laws for the local conservation of mass, linear momentum and kinetic energy 
associated to the Boltzmann equation, one must consider velocity averages of solutions of (1) with 
respect to an appropriate integrable function 0 : R'^ —> R, and in turn use elementary properties of 
the family of Boltzmann scattering matrices {crnlnggrf-i. Indeed, one can show formally that 

^ r <f>fdv + Vjc-f(f>fvdv = ^f Cif,f)i(!>iv) + (l>iv)-fiv'„)-(l)iv„))dv, 

Ot jRrf 4 JrJ 

whence 

^ r ffdv + V^- f ffvdv^O 

ot jRrf jRrf 

if (p satisfies the identity 

fiv'n) + 0(V^) ^ (Piv) + 0(V), (3) 

for all V - [v, v] e R^*^ and n € By choosing f - (p{v) to be 1, v or ^Ivl^, one recovers PDE 
expressing the local conservation of mass, linear momentum and kinetic energy for /, respectively. 
Another important observation in the theory of the Boltzmann equation is that the entropy 


map 


/ 


f f 

Jr‘‘ 


f log / dxdv 


is a formal Eyapunov functional for the dynamics generated by (1), since it can be shown that 


[ 

Jr*' 


C(f, f) log fdv = 


= -- f 

4 Jr 


r c(/,/)iog 



i // J 


dv<0 


( 4 ) 


with equality holding if and only if / is a Maxwellian distribution /m. 


/m(v) ^ 


P 


( 2 n - 0)^/2 


exp - 


| V - u\^ 
20 


for some p, 0 > 0 and u € 


In order to demonstrate that all minimisers of the entropy functional (4) (amongst a suitable 
class of admissible functions) are indeed Maxwellia, one also needs to characterise all solutions 
^ : R^^ —> R of the functional equation (3). Knowledge of all collision invariants also provides us 
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with the nullspace of the linearisation of the collision operator (2) about a global Maxwellian 
/m, which is crucial when it comes to investigating the behaviour of perturbations of equilib¬ 
rium solutions of the Boltzmann equation (1). Moreover, characterisation of collision invariants is 
important for establishing rigorous connections between the Boltzmann kinetic equation and the 
Euler and Navier-Stokes equations of fluid dynamics: see Bardos, Golse and Levermore [3,4] for 
more on such ideas. 

Under various assumptions on 0, it has been shown in the work of many authors (for instance 
Boltzmann [5] for the case, Gronwall [11] for the C° case, Cercignani [6] for the Maxwellian- 
weighted case, and Arkeryd [1] for the case) that if a scalar function 0 : — > R. satisfies 

0(v[,) -I- 0(v^) = (p{v) + 0(v) for all U = [v, v] e and n € it is necessarily of the form 

0(v) = a + b ■ V + c|vp, 

for some constants a, b \,..., bd, c e R.. Any such function (p is known as a collision invariant, as the 
value of the map \y,v\ (p{v) + (p(y) does not change when ‘pre-collisional velocities’ are changed 

to their ‘post-collisional’ values by a-„ for any n e In this article, we will focus our efforts 
on establishing the analogue of this result when the particles in the underlying dynamical system 
are no longer perfectly spherical. 

Although the motivation for studying collision invariants can be found at the kinetic level, we 
make no further study of the Boltzmann equation in the sequel. In all that follows, we focus our 
attention solely at the level of particles. 

1.2. Informal Statements of Main Results. As it takes quite some effort to set up precise 
statements of the main results of this article, we state them at first in a somewhat informal manner. 
For simplicity, we work in two spatial dimensions in all the sequel, i.e. we consider the motion of 
two-dimensional particles evolving in the whole space R^. However, all results in this article can 
be extended to the case of three-dimensional particles evolving in the whole space R.^. 

We study the dynamics of systems of non-spherical particles P consisting of two identical 
compact, strictly-convex subsets with analytic boundaries, i.e. dP is of class C". Naturally, we 
stipulate that at no time should the particles interpenetrate. As such, we must construct a dynamics 
on a suitable phase space of hard particles (see section 3 below for the precise definition of ‘hard 
particle phase space’). The dynamics of the hard particles is governed by Euler’s Eaws of Motion, 
the analogue of Newton’s Eaws for continuum rigid bodies. The first result of this article concerns 
the existence of solutions to Euler’s equations for their evolution which conserve the total linear 
momentum and kinetic energy of initial data for all time, and which also ensure non-penetration 
of the particles for all time. Informally stated, we establish the following result: 

Theorem 1.1. Consider two identical compact, strictly-convex particles with analytic bound¬ 
ary. There exist global-in-time classical solutions to Euler’s equations of motion on the phase 
space of all particle configurations for which there is no particle interpenetration. Moreover, these 
classical solutions conserve the total linear momentum and kinetic energy of any given initial 
datum for all time. 

The precise version of Theorem 1.1 is stated as Theorem 3.1 below. The proof of this result 
makes use of the general existence theory of Ballard [2] for dynamics of rigid bodies with non¬ 
penetration constraints. However, in order to invoke his theory one must first construct scattering 
matrices which resolve collisions between two compact, strictly-convex sets in such a way that total 
linear momentum and kinetic energy are conserved. This construction is performed in section 3.3 
below. The reader might notice that the statement of Theorem 1.1 does not claim that total angular 


4 


L. SAINT-RAYMOND and M. WILKINSON 


momentum is conserved for all time by classical solutions (whose precise definition is given in 3.3 
below). In fact, it has been shown in Wilkinson [17] that classical solutions of Euler’s equations 
which conserve total linear momentum, angular momentum and kinetic energy of initial data for 
all time do not exist for all possible initial data. It is for this reason we confine our affenfion in 
fhis arficle fo dynamics which conserve only linear momenfum and kinefic energy, since fhe notion 
of scaffering map and classical solutions fo Euler’s equations of mofion are infimafely relafed fo 
one anofher. Eef us also draw affenfion fo fhe facf fhaf if may, af firsf glance, seem fhaf our choice 
of dynamics is somewhaf arbifrary, since one can consfrucf disfincf families of solufion operafors 
[rdfgR associated fo Euler’s equations which conserve fofal linear momenfum and kinefic energy 
for all time. We justify our particular choice of dynamics [Td/eR in section 3.5.2 below. 

While fhe spatial collision configuration of fwo hard disks can be characterised by fhe single 
angle fhaf fhe line connecting fheir cenfres of mass makes wifh a given reference line, we nofe fhaf 
an elemenf /3 of fhe fhree-forus is required fo characferise fhe spatial collision configurafion of 
fwo compacf, sfricfly-convex parficles which are nof disks. To see fhis, one mighf wish fo consulf 
figure 2 below. Wifh fhis in mind, we presenf an informal sfafemenf of fhe main resulf of fhis 
article. 

Theorem 1.2. Suppose a measurable map ^:R^xRx§^— satisfies the functional identity 
for collision invariants given by 


(piy'p, Up, d) + (p{v'p. Up, ■&) = ip{v, u, &) + (p{v, u, ■&) 


for every V = [v,v,u,cJ] e R^ and all jd € T^, where [v'p,v'p, u'p,7lfp\ € R^ denotes the post- 
collisional values of the vector V corresponding to the spatial configuration fi. Then ip is neces¬ 
sarily of the form 



for some constants b\,b 2 ,c and some measurable function a : —> R. 

The precise sfafemenf of fhis resulf appears as Theorem 4.1 below. 

1.3. Structure of the Article. In section 2, we revisit the case of hard disk scattering and 
present a new proof of characterisation collision invariants. We derive the equations of motion for 
the physical evolution of hard particles in section 3. The concept of scattering map and regularity 
of solutions of Euler’s equations are intimately linked, so in sections 3.3 and 3.5 we construct 
families of scattering maps and, in turn, classical solutions to Euler’s equations of motion. In the 
final part of the paper, namely section 4, we characterise collision invariants for compact, strictly- 
convex non-spherical particles. New results by C. Viterbo on generators of orthogonal groups of 
matrices, which allow us to establish the proof of Theorem 1.2, are stored in the appendix A. 

2. Characterisation of Collision Invariants for Hard Disks: A New and Simple Method 

Before we embark upon the problem of characterising collision invariants for general convex 
particle scattering maps, it will be helpful to recall the theory which has been established in the 
case of spherical particles (or, more appropriately in our two-dimensional setting, particles which 
are disks). Our approach to this problem appears to be new, and has the advantage of requiring 
no regularity or integrability conditions on the collision invariant <p, only that it be measurable. 
Although we only discuss scattering of hard disks in R^ in this section, all our results also hold for 
the scattering of hard spheres in R^. 
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2.1. State-of-the-art of Previously-established Results. For any S \ consider the asso¬ 
ciated Boltzmann scattering map cr^ : R"* —> R^ for two hard disks given by 


(T^VV] := (/ - 2r^ V with F = [V, F] e R^ 


(5) 


where 

e(<A) 

-e(iA) 


_ 1 
7^ —p 

V2 


with e{4f) := (cos sin i/^) € R^ and tjj denotes the angle that the line connecting the centres of 
mass of the colliding disks makes with the positive v-axis. One can check that for every choice of 
(/^ € the scattering map cr^j, conserves total linear momentum, angular momentum and kinetic 
energy of any given velocity vector V € R^. Under the assumptions that 0 : R^ —> R be in Lj'^^(R^) 
and satisfy the functional equation 


+ 0(v^) = 0(v) -I- 0(v) (6) 

pointwise almost everywhere on R"* x S', where the post-collisional velocities and are given 
in terms of cr^[F] = (a-^[F]i, ...,cr^[F] 4 ) as 


V 




CT^[V]l \ 

I 


and 


I ’ 


it has been shown by Arkeryd ([1], lemma 2.8) that (p is necessarily of the form (p{v) = a+b-v+c\v'^ 
almost everywhere for some constants a,b\,b 2 ,c e R^. Our new proof of characterisation of 
collision invariants covers the case where (p is only measurable on R^, as opposed to being of class 
Lj'^^(R^). On the other hand, we ask that the identity (6) hold for all i/^ € S' and for all V e R^. In 
order to produce the most general result possible, one would need to extend our argument to the 
case where (6) holds for almost every i/^ € S' and almost every V e R^, as opposed to everywhere 
on s' and R"^, respectively. We do not attempt do this here. 


2.2. Orbits of Scattering Groups on R"'. In order to motivate our new group-theoretic ap¬ 
proach in the case of general strictly-convex particles, let us rewrite identity (6) as 


= %{y) 


(V) 


for F € R^' and ijj e S', where 

^ci>(y) := (p{v) + (p(y), 

with F = [v, v] € R'', assuming that cp be only measurable and, without loss of generality, that 
0(0) = 0 and thus <1>0(O) = 0. In particular, identity (7) implies that for any fixed choice of F and 
any collection of angles 0i,..., 0^. € S', one has 

° ■■■ ° 

namely that the map is constant on the left group orbits GF c R"' for any given F e R"', where 
G c 0(4) is the group generated by the 1-parameter family of reflection matrices {/ - 2y^ ® : 

0 € s'}, namely 

G:-({/-2y(,®r^ : 0eS')). (8) 

Let us now find the group orbits GF for any F e R^. For e > 0 and p e R^ satisfying > |pP/2, 
we define M(e, p) to be the subset of R"^ given by 


/ 7i + Fa 

\ Fa + F4 



M(e, p) ] F e R^ : |Fp - and 
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which is evidently homeomorphic to When = |pP/2, M(e, p) is a singleton and when 
< lpP/2, one can check M(e, p) is empty. It is clear that when V e R.^ is given, the Boltzmann 
scattering matrix cr^ maps M(e, p) to itself for any i/^ € where e = | V| and p = (Vi + Vj„V 2 + V 4 ). 


2.3. Reduction to Canonical Form. As the sets M(e, p) are homeomorphic to for > 
IpP, we can expect to reduce our study of scattering groups acting on M(e, p) to the study of some 
other group acting on To show this, we reduce our problem to a kind of canonical form. Indeed, 
for > lpP/2, we consider the bijection /ie,p : M(e, p) —> given by 


with inverse given by 


Fi- 1/3 

V(^l - ^3)2 + (^2 - F4)2 I ^2 - F 4 


for V e M(e, p). 


- 2 


V2e2^]^^i + Pi 

V2eMpR^P2 

Pi - 

P 2 - V2e2^n^^2 


for ^-(^i,^ 2 )eSP 


One has that cr^ e maps F to (/ - 2y^ ® if and only if the matrix 


■= I - 2e{il/) ® e{i}/) € 


5)2x2 


maps /ie,p[F] to (/-2e(i/^)®e(i/^))/ie,p[F]. Thus, if the group ({ 5 ^ : ip e SM) £ 0(2) acts transitively 
on the circle Sp it will follow immediately that the group orbit GV is identically equal to M(e, p). 
This is indeed the case, as the following elementary result shows. 


Proposition 2.1. The group {{I — 2e{ip) ® eip) : ip & S^}) Q 0(2) acts transitively on 


Proof. For any two points ^1 = e{ip\) and ^2 = e{\p 2 ) for ipi,ip 2 ^ we set ip' {ip\ + 
ip 2 )l'^ £ s'- One can check that C 2 - {1 - 2e(iA')"'" ® ^(<AO"'")^i> and so we are done. □ 

Transforming back to 1.^, we immediately infer that the orbits of points of under the action 
of the scattering group G in ( 8 ) above are given by 

( M(e,p) if >^- 

GV^l 

\ j[ 2 p, 2 p]) if e 2 = ^ 

Since is constant on each left orbit GV, it follows that 

%{V)^%{v + v, |vp + |vp) 

for some new measurable function < 1)0 : R^ x R —» R. One may then check (using the fact that 
0 (0) = 0) that < 1)0 satisfies the identity 

I’.aCa |vp) + <l) 0 (v, |vp) ^ (l) 0 (v + V, |vp + |vp) for all F = [v,v] € R"^. (9) 

It is at this point we appeal to results on the characterisation of solutions to Cauchy’s Functional 
Equation (see, for instance, the book of Kuczma [13]). 


2.4. Results on Cauchy’s Functional Equation. We recall that, under the assumption / : 
R^ —> R be a measurable function, any solution of the functional identity 

fix) + fly) = fix + y) for all x, y e R^ ( 10 ) 

is necessarily of the form fix) = cx for some c e R. We remark in passing that one cannot weaken 
the assumption that <p is measurable, if one wishes to avoid dealing with ‘pathological’ solutions of 
Cauchy’s functional equation. Indeed, by dropping the assumption of measurability and assuming 
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the axiom of choice, it has been shown by Hamel [12] that there exist discontinuous solutions of 

(10). 

One can use the fact that all measurable solutions of (10) are of the form f{x) - cx to 
characterise all measurable maps satisfying the functional equation (9) for above. We now 
quote a result contained in Truesdell and Muncaster ([15], pages 72-73 and pages 88-89), whose 
proof we revisit in detail in section 4.5. 

Proposition 2.2. Suppose that a measurable map <1) : K.^ x K. —> R satisfies the identity 

<1>(V, |v|2) + (1)(V, |v|2) = (I)(v + V, |vp + |v|2) (11) 

for all V, V e R^. It follows that d) is necessarily of the form d)(v, |vp) - b ■ v + c|vp for some 
constants b\,b 2 ,c 

Using the additional observation that any constant function is also a collision invariant, it 
quickly follows that if a measurable function 0 : R^ —> R satisfies the identity 

0(vp + - (p{v) + 0(v) for all v € R^ and i/^ € 

then it is necessarily of the form 0(v) = a + b ■ v + c|vp. As such, one can view the problem 
of characterisation of collision invariants as the problem of classifying all scalar invariants of a 
given group action (namely that of the scattering group G) on Euclidean space R"^. It appears that 
this perspective on the problem is new. In particular, we emphasise that we placed only minimal 
assumptions on f, namely that it be only measurable on R^. It is this group-theoretic perspective 
on the problem we adopt in order to prove the main result of this article, namely Theorem 1.2 
(restated precisely as Theorem 4.1 below). We now leave the case of hard disks to study general 
compact, strictly-convex sets with C" boundaries. 


3. Dynamics of Compact, Strictly-convex Particles 

Although collision invariants themselves have no relationship to particle dynamics, what con¬ 
stitute pre- and post-collisional velocities at collision is, however, inherently a dynamical issue. It 
is for this reason we must address the dynamics of particles in this article. As collision invariants 
only involve two-particle interactions, we study in all the sequel the evolution of two compact, 
strictly convex sets t i-^ P(t) and t i-^ P(t) in the plane R^ which do not interpenetrate. We as¬ 
sume that their boundary curves are of class C", and that the motion of P and P takes place in 
the absence of external forces. We subsequently refer to compact, strictly-convex subsets of R^ as 
hard particles. As there are no externally-imposed forces in our systems under consideration, the 
evolution of the sets P(t) and P(t) before collision is determined by their initial states, namely their 
initial spatial configurations (centres of mass and orientations) and initial velocities (both linear 
and angular). In order to construct a ‘physical’ evolution for these two hard particles on R^, we 
appeal to Euler’s Laws of Motion for continuum rigid body classical mechanics. We recall that 
Euler’s laws are the appropriate extension of Newton’s laws of motion to the study of continuum 
rigid bodies. We refer the reader to Truesdell ([16]) for more on this topic. 

Eet us now set up the basic objects with which we work throughout this article. Suppose that 
P* c R^ is a compact, strictly-convex set with boundary of class Moreover, suppose that its 
centre of mass lies at the origin, i.e. 


£ ydy = 0. 
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We shall subsequently eall any such set a reference particle. When an arbitrary centre of mass 
A e and orientation ?? e S' have been given, we write the v-translate and ?9'-rotation of P* as 

P(a, i9') /?(7?)P* + X, 


where R{a) € SO(2) is the rotation matrix 

R{a) := 


cos a - sm a 
sin a cos a 


The evolution of the sets P(0 and P(t) is expressed by 

P(0 /?W0)P* + x{t) and P(t) /?(^(0)P* + x(0, 

with the centres of mass x{t),x(t) e R^ and orientations i9'(t), &{t) e S' being related to the linear 
velocities v(t),v(t) € R^ and angular speeds € Rby the formal differential relations 


together with 


dx , dx _ 

— =v and “T - 

dt dt 


dd dd _ 

-— = o) and -—= CO. 
dt dt 


( 12 ) 


(13) 


We gather the spatial and velocity data into single phase vectors z and z given by 
z(0 - Vxit), d{t), vit), co(t)] € A1 R^ X s' X R^ X R, 


and also 

Z(0 - [x(t), d{t),v(t),co{t)] € At R^ X s' X R^ X R. 

We define fhe single phase vector which characferises fhe sfafe of fhe whole sysfem af fime t e R by 
Z(t) := [z(0>z(0] £ As we do nof wish fhaf P(t) n P(t) have posifive 2-dimensional Lebesgue 
measure for any fime t, we sfipulafe fhaf fhe range of fhe maps t Z{t) belong to fhe phase space 
2)2 = d) 2 {P*) defined by 

n 2 iP*) := jz € : card P{x,d) n P{x,d) < l), 

where Z = [z, wifh z = [x, d, v, io\ and z - [a, d, v, iZi]. As if will be useful in whaf follows, we 
define fhe associafed spatial projecfion operafor fli : D 2 —> R"' x by fhe rule 

fliZ [x,l,,d,d] when Z = [z,z] e D 2 . 

We also define fhe velocify projecfion operafor 112 : ID 2 —> by fhe rule 

n 2 Z [v,v,co,cJ] when Z = [z,^ € D 2 . 


In order to be complefely correcf, we nofe fhaf fhe differenfial relations (12) and (13) only hold 
in general af fhose fimes t € R for which P(t) n P(t) = 0, i.e. fhe fwo-sided derivative limifs in 
(12) and (13) hold af fhose fimes t when P(t) and P(t) are nof in collision wifh one anofher. Af fhis 
poinf, if will prove helpful to make fhe following definition. 

Definition 3.1. For any Zq € ID2 and a map Z : R ^ ID2 safisfying Z(0) = Zq, we define fhe 
associafed sef of collision times T (Zq) c R to be 

T{Zf) := [t e R : card P(0 n P(0 = l). 


9 


In order to derive the equations of motion which govern the particles P(t) and P(f), we first 
of all consider a class of 1 -parameter families of operators {r,},gR (r, : 1 D 2 —> 2^2 for each t € K.) 
for which the maps t 1 -^ YViTtZ^ and t i-> n2r,Zo have ‘reasonable’ analytical properties. Indeed, 
in order to make concrete the primary objects of interest in this article, we make the following 
important definition. 

Definition 3.2. We shall call a family of operators {r,),gR with Tt : D 2 —> 2)2 for each t € R 
a hard particle flow on D2 if and only if for any Zq e D2, the map t 1-^ nir^Zo continuous and 
both left- and right-differentiable on R and the map t i-> n 2 rfZo is lower semi-continuous and left- 
differentiable on R. Moreover, we stipulate that both 1 1 -^ Ilir^Zo and 1 1 -^ n 2 r(Zo be differentiable 
at all times t for which TfZo € D 2 \ 52)2- 

The class of hard particle flows on D 2 is evidently a rather large one. A basic question in 
classical mechanics is the following: “Which hard particle flows on D 2 can one consider to be 
physical!” To answer this question, and to specify in precise mathematical terms what we mean 
by physical, we appeal to Euler’s Laws of Motion. When deriving an appropriate set of DDEs 
that govern the evolution of the phase map t i-> Z(t), we divide our considerations into two cases, 
namely those times during which the dynamics is collision free, and those times at which a collision 
takes place. 

3.1. Deriving the Equations of Motion when P(?) n Pit) - 0. Suppose a hard particle 
flow {EjlfeR on the phase space D 2 has been given. This flow gives rise naturally to a map U : 
R^ X R X 2)2 —> which provides the instantaneous linear velocity of any material point x in 
R^ at any time t, once an initial condition Zq e D2 has been provided. Indeed, recall that if the 
centre of mass x(t) of a planar rigid body P(t) translates with linear velocity v(t), and P(t) rotates 
with angular speed cuff), then the linear velocity of any other point on the body is expressed by the 
formula 

v(y, t) = v(0 -I- ojiOiy - x(t))-^ for y e R{&{t))P^, + x{t), 

where y-*- := (-y 2 ,yi) for any given y = (yi,y 2 ) e R^- As such, the map U is given explicitly in 
terms of {r,|fgR by 


v(t) -I- - x(t))-‘- if X € P(t), 


Uix,f,Zo) - ■ 


v(0 -I- (joit){x - x(t))-‘- 


if X e P(0, 


0 


otherwise. 


where Ifir^Zo = [x(0, x(t), i2(t), and n 2 rfZo = [v(t), v(i),m(0,m(0]. Since {Tt}teR is a hard 
particle flow (definition 3.2), it follows that t U{x, t,Zo) is a differentiable function at all t for 
which TiZq € 2)2 \ (92)2- 

We appeal to Euler’s laws of motion in order to partition the class of hard particle flows into 
‘physical’ and ‘unphysicaT flows. We henceforth assume that the motion of the hard particles P 
and P takes place in the absence of external forces. Consider any Zq € D 2 for which 7~ (Zq) R, 
and let us restrict our attention to the open set /(Zq) c R on which P(t) n P(t) - 0, i.e. where 
the map 1 1 -^ 1122 ,Zq is differentiable. We now consider Euler’s First Law of Motion (Truesdell 
[16]), which states that for any smooth evolution of smooth subsets 1 1 -^ Q(t) c R^, a physical hard 
particle flow should satisfy 

± f 

dt 


Uix,f,Zo)dx - 0 for t e l{Zo). 


(14) 
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Since we are free to choose the family of testing sets {Q(t) : t € /(Zq)) as we wish, we first pick 
it to be a family of smooth open sets such that P(t) c Q.{t) together with Q(t) n P(t) = 0 for 
all t € /(Zo). Similarly, we can also choose Q.{t) to contain particle P(t) alone. As t i-> TtZ^ is 
differentiable on /(Zq), identity (14) reduces under these two choices to the ODEs 


m — = 0 and m — = 0, (15) 

dt dt 

where m - Jp dy is the mass of the reference particle P*. This implies in particular that the total 
linear momentum of the initial datum Zq is conserved on /(Zq). Thus, in the absence of external 
forces and collisions, Euler’s first law simply reduces to the conservation of linear momentum. 

It is now we turn to Euler’s Second Law of Motion (Truesdell [16]), which states that 


d r 

— (x- a) ■ U(x, f, Zq) dx = 0. 

dt Jn(i) 

By appropriate choices of Q(t), we discover that Euler’s second law of motion reduces to 


—( - m{a - x(0)"*" • v{t) + J(o{t)) = —{ - m{a - x(t))'‘~ ■ v{t) + Joj{t)) = 0, 
dt dt 


where J Jp \y^ dy is the moment of inertia of the reference particle P*. By appealing to the 
ODEs (15) derived above, we may infer that 



and 



Therefore, it is clear that Euler’s first and second laws together imply the conservation of linear 
and angular momentum for T^Zq on /(Zq). Importantly, one may check that Euler’s first and second 
law imply that total kinetic energy is conserved in time, in the sense that 


£ f 

dt Jr2 


\U{x, t;Zo)P dx -0 


for all t € /(Zq). 


With this discussion in place, we now specify in precise terms what we mean by a classical solution 
to the ODEs derived from Euler’s laws. 


3.2. A Dynamical System and its Boundary Conditions. Due to the possibility of particle 
collisions, we cannot expect the velocity maps t Yi^fTtZo to be differentiable both on the left 
and on the right on R. As such, we separate out the information contained in Euler’s ODEs into its 
left- and right-limits. We consider the following class of dynamical system, namely the evolution 
of two identical compact, strictly-convex sets P and P (which are translations and rotations of the 
reference particle P*), whose phase trajectory t i-^ Z{t) e ID 2 satisfies the system of one-sided 


ODEs 


d_ 

dt- 


in the classical sense for all t € R, where 

x{t + h) - x(t) 


X 


V - 


X 


V - 

d 


tx)- 

. d 

d 


oj - 


— 


and — 


— 


V 


0 

1 

V 


0 

tx) 


0 


u 


0 


v-(t) := lim 
h-~^0- 


h 


and co-(t) := lim 
/ 2 —> 0 - 


■&{t + h) - df) 
h ’ 


and similarly for the barred variables v_ and m_. We also ask that t Z{t) satisfies the system 


X 


v+ 


X 


V+ 

d 


Cx>+ 

and — 

d 



V 


0 

dt+ 

V 


0 

tx> 


0 


oj 


0 
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in the classical sense for t e 1. \ T"(Zq), where 

x{t + h)-x{t) 
v+m := lim — 

/i^0+ 


and 


— w +(0 := lim 

h /i^o+ h 

and similarly for the barred variables. With this in place, we make the following definition. 


d-{t -\-h) - {)-{t) 


Definition 3.3. For a given initial datum Zq e D 2 , we say that Z : R ^ D 2 is a (global-in¬ 
time) classical solution of the above system of Euler’s equations of motion if and only if v, v, •& 
are continuous piecewise linear on R, and v, v, a>, to are lower-semicontinuous piecewise constant. 
Moreover, these maps satisfy the above DDEs pointwise on R for the left-derivatives, and pointwise 
on R \ 7” (Zq) for the right-derivatives. Einally, Z(0) = Zq. 


Evidently, the system of ODEs above is not enough to determine a family of flow operators 
{r,)fgR on 2)2 uniquely. Indeed, one must specify how to update the dynamics at all collision times 
T € 7”(Zo), i.e. for all t such that 

card P(v(t), 77(r)) n P(x(t), &(t)) = 1. 


When two compact, strictly-convex nonspherical particles are in contact at a single point, their 
configuration can be characterised (with respect to the reference particle P*) by an element /? of the 
3-torus T^. In order to be able to construct a flow on 2)2 ^ one must in turn construct an associated 
family of velocity scattering maps on R®, each member of which sends elements in a set 

of ‘pre-collisional’ velocity vectors to elements in a set of ‘post-collisional’ velocity vectors (see 
section 3.3.3 below for the precise definition of these sets). 

Not only this, one would ideally wish the family of flow operators {TfjfgR on D 2 to conserve 
the total linear momentum, angular momentum and kinetic energy of any given initial datum; 
consequently, any scattering map oy? : R^ —> R® should also have this property (consult section 
3.4 for a discussion of the conserved quantities of the dynamics). However, it is shown in [17] that 
such a family of scattering maps on R® does not exist. Nevertheless, we study a class of scattering 
matrices which gives rise to a hard particle flow on 2)2 that conserves total linear momentum and 
kinetic energy of all initial data. Before we can construct any flow associated with the above ODEs 
on 2)2 ^ we must first find a convenient way by which to parameterise collision configurations. This 
is the subject of section 3.3.1 below. 

3.3. Scattering Maps on R®. Scattering maps are the fundamental objects with which we 
work in this article. In particular, they must be constructed if one is to employ the existence theory 
for rigid body mechanics due to Ballard (see, in particular, hypothesis H3 [2] p.212). In order to 
construct scattering maps, we must first find a careful parameterisation of all possible two-particle 
collision configurations, and then in turn specify what one means by pre- and post-collisional 
velocity vectors. 

3.3.1. Parameterising Collision Configurations. We now parameterise the set of all Z e D 2 , 
up to translation, such that card P{x, ■&) n P(x, i?) = 1. In this direction, we consider what we call 
a reference collision configuration which will allow us to parameterise a general collision config¬ 
uration of two particles by an element of the 3-torus T^. By considering the plane R^ furnished 
with polar co-ordinates, we make the problem of describing collision configurations considerably 
simpler. Indeed, as previously indicated, it will be of some help to consider the centre of mass of 
the reference particle P* as at the origin of R^, which the polar map 

i ip cos when (p, i/^) e (0, 00 ) X 

( 0 , 0 ) otherwise. 
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co-ordinatises. We shall use P* to define reference collision maps which are functions of the 
polar angle iff € and the orientation 0 € S' of the particle exterior to the reference particle 
P*, namely ne = ne{tfr),Ng = Ne{tJ/),pe = P 0 {>Jj),qe = qeW and dg = dg{iJ/)\ see Figure 1 below 
for an illustration of these quantities. They constitute the essential spatial data used to construct 
post-collisional velocities in a collision between two particles. 

We begin by making the following definition. 

Definition 3.4. Let 0 € S' be given. The distance of closest approach dg{ij/) of the centres 
of mass of P* and P(-, d) for the given elevation angle \j/ is defined to be 

dg{if/) ■.= inf |(i > 0 : card P* n (7?(0)P* + de{4/)) = o|, 

where eiijj) (cos ifj, sin ifj). 

With this basic and important quantity defined, we make another important definition. 

Definition 3.5. We say that two particles Pi, P 2 c are in a reference collision config¬ 
uration whenever P,- = P* (for some i € {1,2}) and there exist 0 € S' and 1 /^ € S' such that 
P j = /?(0)P* + dg{f)e(f/), for j 4 i. 

The other basic collision configuration quantities are now straightforward to characterise. We 
define the collision vector p = pe(<A) to be the unique element of the set 

P* n (/?(0)P* + dg{f)e{f)), 

and the conjugate collision vector q = qgif) by 

qgii//) := pg{ilf) - dg{f)e{f). 

Since 5P* = P* \ int(P*) is a closed C" curve in R^ and can therefore be described locally by a 
smooth polar map /ip., one can speak of the (outward) contact normal n = nglfj) to the point of 
collision p = pgiijj), which is given by 

h'piagif))^ 

■ |/ip^ (o'e(iA))-^r 

where 

ag{f) — arctanf^^^'j. 

The exclusion normal Ng = Ngifj) is defined to be the (outward) unit normal to the closed C" 
curve Cg given by 

Cg -.= {<ie(iAMiA) : lA £ s'}. 

Notice that in the case of hard disks (when P* = B(0, r) for some r > 0), this curve is simply a 
circle of radius 2r, whence Ng coincides identically with ng. These basic vectors are illustrated in 
Figure 1 below. 

3.3.2. General Collision Configurations. When two particles P and P in the dynamical sys¬ 
tem described above satisfy card P(t) n P(t) = 1 for some r € R, we shall say they are in a general 
collisional configuration. Of course, it is not the case that they are necessarily in a reference colli¬ 
sion configuration as defined above in definition 3.5. In order to solve for the post-collisional linear 
velocities and angular speeds of two particles with arbitrary orientations (described by i9', i? e S') 
and arbitrary relative position (described by (A e S'), it is expedient to relate general collisional 
configurations to the reference configuration introduced above. 
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Figure 1. An example of a reference configuration for P* and P = /?(0)P* + 

If P* remains the standard reference particle, suppose P, P are of the form 

P = /?(j?)P* and P = 7 ?(j 9')P* + x, 

with X e K.^ such that card P n P = 1, i.e. P and P are in a collisional configuration. Thus, there 
exists an angle of elevation i/^ € S ^ and a constant g - Q{d, > 0 such that 

P = R{d)P* and P -/?(^)P* + (16) 

In order to write down the appropriate distance of closest approach d ^, together with the analogous 
collision vector p ^, its conjugate and the normals and in terms of the respective quantities 
dg , pg , qg , ng and Ng defined above, we perform some rotations. Acting on the system described in 
(16) by the rotation matrix 

cos ■& sin ■& 

- sin ■& cos d 

we map P to P* and P to Rid - j?)P* + g{d, d, - d). This transformed system is now in a 
reference collision configuration. In particular, Q{d, d, ijf) = - d). Finally, by rotating back 

to the original configuration described by (16), it is clear that the basic collision quantities for two 
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identical particles of orientations ■&,& € whose centres of mass define a line of elevation tp with 
respect to the polar axis are the following: 

dp = d^{ip) d-^_^{>p - ■&) (distance between centres of mass) 

Pp ^ ■= R{'9)p:ff_^{'P - d) (collision vector) 

qp = R{d)q-^_^,^{>p - d) (conjugate collision vector) 

Np = N'liP) R{'d)N-^_^{ip - d) (exclusion normal) 

and 

np = n^((A) := R{d)n:^_^{ip - d). (outward contact normal) 

These are illustrated in Figure 2 below. We work with these five fundamenfal vecfors in all fhe 
sequel. 

Remark 3.1. As we have done above, we shall offen wrife fhe quanfifies such as d'^{ip) simply 
as dp wifh [3 = {d, d, tp) when fhe values of i9', ??, e S' are undersfood. If will offen be convenienf 
fo use fhe nofafion d^{p) whenever we emphasise fhaf fhe paramefers {d, d) e have been fixed, 
and p i-> d'^ip) is considered a funcfion of p alone. In fhis case, when fhe values of {d, d) € are 
undersfood, we shall simply wrife d{p). This allows us fo make fhe presenfafion of our argumenfs 
(especially fhose in section 3.6) less cumbersome. 

3.3.3. Pre- and Post-collisional Velocities in R^. We now consfrucf scattering maps on R^ 
which assign posf-collisional velocifies fo pre-collisional velocifies of fwo parficles in a collision 
configuration in such a way fhaf 

card P(x(f), d{t)) n P(x{t), d{t)) < 1 

for all t in a sufficienfly-small neighbourhood of a given collision fime t € 7~ (Zq). Once we have 
such a map fhaf uniquely updafes fhe particle velocifies, we may consfrucf a global flow on phase 
space Dj corresponding fo a classical solufion of fhe sysfem of governing ODEs infroduced in 
secfion 3.2 above using fechniques from [2]. 

We now derive sefs of pre- and po^t-collisional velocify vecfors, and define whaf we mean by 
a scaffering map. In order fo do fhis, lef us consider fhe auxiliary map F : R'' x ^ R given by 

F{x, X , d, d) := \x-x\- cTJ| arcfan 

Clearly, F{x,x, d,d) > 0 if and only if P(x, i9')nP(x, d) = 0 ; moreover, F{x,x, d,d) = 0 if and only 
if card P(x, ??)n P(]c, d) = 1. We now infroduce a hard particle dynamics {TflfgR associated wifh fhe 
ODE sysfem in secfion 3.2 above. Consider fhe maps {x,d) : R —> R^ x S' and (x, ?9') : R —> R^ x S' 
(wifh niTfZo = [x, X , d, ??]) which satisfy 

card P(x(f), d{t)) n P(x(0, d{t)) < 1 

for all time t e R. We recall fhaf, by assumpfion, Z(t) = IIi TfZo is bofh leff- and righf-differentiable 
af all times, fhe only poinfs af which righf-derivafives do nof necessarily agree wifh fhose on fhe 
leff being fhe sef of collision fimes 7~ (Zq). 

Consider now any collision fime r € 7” (Zq). Using fhe assumpfion of leff-differentiabilify of 
fhe relevanf phase maps, we have 


X2 - X2 

Xi - Xi 


d 
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Figure 2. A general eollision eonfiguration 


for arbitrary t €T (Zq), whieh a ealeulation reveals to be 

e(iA) - - ~ • V- - e(iA) - - ■ v- 

with and i9'(t) simply denoted by x,x,& and j9', respeetively. Moreover, sinee 

the partieles are in a eollision eonfiguration, there exists ^ satisfying the identity x -x - 
Now, we notiee that the eurve of elosest approaeh 

has (non-normalised) normal veetors 

1 dd-Z a 

:= -Tr^(*A - , 


( 17 ) 
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whose normalisation we denote by := N^(tJ/)/\N^{tJ/)\. Moreover, we make the observation 

that 

We therefore write the above inequality in the more compact form 

ivfW-v.-ivfW-v. 

+ (*A) - (<A)) ^ 0, 

where is the vector 

riifp) - ~ 

As it is one of the most important quantities in all that follows, we make the following definition. 


Definition 3.6. For any yS e T^, the collision normal yp e R° is defined to be 


yp ■= 




Np 

-Np 

[rp - dpeitl/)Y ■ Np 




■p ^yp 

where 

2 1 I -L 

• Np 

Remark 3.2. A quick calculation reveals that the collision normal yp is not of unit norm. It 
will be useful rather often to employ the unit collision nomial'yp M~^yp in what follows. 


+yp 


■Np\\ 


(18) 


In the language of definition 3.6, one then has that 


if and only if 


—F{x(t), x(t), d{t), d{t)) 
dt- 


yp-V- < 0, 


< 0 


where V- - [v_, v-, m-, w-]. In a similar way, one can treat the post-collisional case and deduce 
that 

d _ - 

—F{x(t), x{t), dit), d{t)) > 0 

dt^ t=r 

if and only if 


yp-V+> 0. 

Let a spatial configuration point /? € be given and fixed. Wifh fhe above discussion in mind, we 
define fhe sef of pre-collisional velocifies associafed wifh fhe spatial configuration /3 € fo be 


: F-r^<0), 

and fhe sef of all post-collisional velocities fo be 

: F-y/j >0). 

Evidenfly, U We denofe fhe infersecfion Z^ D Z^ of fhese fwo half-spaces by Z°. Wifh 

these definitions in place, we can now say what we mean by a scattering map on R®. 


Definition 3. 7. We say that a bijective map cr^ : R® —> R® is a scattering map corresponding 
to the spatial configuration /? e if and only if (Tp(L^) - Z^ and o-p o crp = l on M.^. 
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and recalling that jp = M ^jp, we can recast the above 
Pp{P^■%>0, ( 21 ) 

( 22 ) 


Suppose p € T^, i.e. let the orientations and centres of mass of two particles in a collision 
configuration be given, and let crp be an associated scattering map. By definition, 

V-yp<0 ^ (Tp[V]-yp>0, (19) 

and also 

V-yp>0 ^ (Tp{V]-yp<0. (20) 

It will be convenient in the rest of this article to write the above inequalities in what we shall call 
quasi-momentum variables. Consider the mass-inertia matrix M € given by 

M := diag( V7, V7). 

Writing P := MV for a given V € 
conditions as 

P-yp<0 
and 

P-yp>0 pp[P] -yp^O, 

where the transformed scattering map pp is given by 

ppVP] := M(Tpm~'Pl 

We write the associated transformed set of pre-collisional velocities as and the post-collisional 
velocities as 

There are many involutions cr^ : R® —> R® which map the lower half-space to the upper 
half-space We now specify some conservation laws from classical mechanics, attributed to 
Euler’s laws of motion, which should be respected by the hard particle flow {r,}(6R on phase space 
ID 2 . In parficular, in view of fhe resulfs in [17], we sfipulafe fhaf fhe flow should conserve only 
fofal linear momenfum and kinetic energy of given inifial dafa Zq € D 2 . 

3.4. Derivation of the Algebraic Constraints. Suppose the particles in collisional contact 
P := /?(7?)P* and P := 7?(??)P* -td^(i/r)e((/^) are given, together with their respective linear velocities 
and angular speeds V e with p - {&, ij/). We seek post-collisional linear velocities and 

angular speeds V e such that there is conservation of total linear momentum and there is 
no loss of kinetic energy following collision. In what follows, unprimed quantities will denote 
pre-collisional ones, while those which are primed denote post-collisional ones. 

Adhering to Euler’s first law of motion, we stipulate that the values of the pre- and post- 
collisional velocities should satisfy the conservation of linear momentum, i.e. 


I Vg(y, T)dy+ I v'g(y, T)dy ^ ( v(y, T)dy+ ( v(y, r) dy, 
JP(z(t)) JPair)) JP(z(t)) JP(I(t)) 


(COEM) 


^P(7(t)) 


P(z(r)) 


'P(z(T)) 


which since v(y, t) = v(t) - 1 - a){t)(y - x{t))^ and v(y, t) = v(t) -P a}{t){y - x{t))^ (and similarly for the 
primed variables) reduces to 


mvp -P mvp = mv -I- mv. 


(23) 


We also require that total kinetic energy be unchanged after the collision of the two particles. The 

conservation of kinetic energy takes the form 

r Wp(y, t)I^ dy+f \vp(y, t)P dy=f |v(y, T)f dy-P f |v(y, t)P dy, (COKE) 
Jp(z(r)) Jp(7(t)) Jp{z(t)) Jp(7(r)) 


^P(z(r)) 

which reduces to 


Mv'p?' + Jloj'pf' + = m|v|^ + Jo? + m\vf' + Jco^. 


(24) 
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Figure 3. A locus of closest approach with the exclusion and contact normals 

Expressing the above conservation laws in scattering map notation, we find that (23) takes the form 

/ crpmi + cTpWh \ + ^^3 \ 

\ crpWh + a-/j[F]4 ) \ F2 + V4 / ’ 

while (24) takes the form 

mcrpiVf = \MV\^, 

where V = [v,v, m,m]. As claimed above, in order to prove Theorem 1.1 (or, rather, the more 
precise statement 3.1), we must first construct a family of scattering maps on R^, each 

member of which conserves total linear momentum and kinetic energy. This is the aim of the 
following section. 

3.5. Construction of a Dynamics for Euler’s Equations on D 2 . We now aim to prove the 
following more precisely-stated form Theorem 1.1. 

Theorem 3.1. Suppose P* c R^ is compact and strictly-convex with boundary 5P* of class 
C‘^. For each Zq € 2 )2(P*). there exists a global-in-time classical solution Z{t) = TtZo of Euler’s 
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equations with the property that 

/ (n,T,z„). . (n.rA), \ ^ / (n,z.), + (n,z„)3 \ 

\ (n2T,Zo)2 + (n2r,Zo)4 / 1 (n2Zi,)2 + (n2Zo)4 1 

and 

\MU 2 TtZ 0 f = \MU 2 Z 0 f for all t € R. 

Notice that the above theorem makes no claim on uniqueness of solutions. They are, how¬ 
ever, unique with respect to a fixed family of scattering matrices {cr^l^gTS. In other words, once a 
family of scattering matrices has been chosen and fixed, the classical solutions of Euler’s equations 
constructed using the theory of [2] are unique. As such, we must make a choice regarding with 
which family of scattering maps we wish to work. Since the study of linear scattering maps and 

their corresponding collision invariants is made possible by means of group theoretic arguments 

for subgroups of the orthogonal group 0(6) (see section 4 below), we subsequently focus on the 
case where scattering maps oy? : R.^ ^ K.® are matrices. One could construct solutions of the 
ODEs in the case when the scattering family is a collection of nonlinear maps on R®. We 

do not, however, pursue this idea any further here. 

3.5.1. The case of linear scattering cr^ : R^ —> R®. We establish the following preliminary 
result. 


Proposition 3.2. For a given fi € T^, let crp be a linear scattering map which conserves kinetic 
energy and linear momentum, i.e. (Tp\y^ satisfies (23) and (24) and for all V € R^. Then o'/j is 
necessarily of the form 


(Tjj = M 


-1 


Ex®Ei+E 2®E2 + ^ Ai(fi)Ei(fi) ® Eiifi) -yp®yp 


i=3 


M, 


where E\ = (-^,0,-^,0,0,0) E2 - (0> "^>05 0). is any orthonormal basis for 

span{£'i,£'22^1'''. d,(yS) € {-1,1) andfp is the unit collision normal (3.6). 


Proof It will be convenient to consider the problem cast in quasi-momentum variables as 
introduced above in section 3.3.3. Indeed, given the scattering map crp we define the map pp{P] 

for P e R^. Since crp is linear if and only if pp is linear, we may suppose that 
PpVP^ = PpP for some Rp e GE(6). Moreover, we also infer that pp is an involution on R^, whence 
7?^ = 7. It will now prove useful to consider the spectral structure of Rp. 

We first note that since the conservation of kinetic energy (24) implies that \RpP\^ - |Pp for 
all P € R^, it follows that Rp € 0(6). Moreover, Rp can only have real eigenvalues A with |/1| = 1. 
Now, the conservation of linear momentum 


/ crp[V], + (TpWh \ + f's \ 

\ (Tpm 2 + crpm^ I { V 2 + V 4 j 

implies that 

RpP • El - P • El and RpP • E2 - P • E2 for all P € R^ 

where Ei = (1,0,1,0,0,0) and E2 = (0,1,0,1,0,0). We immediately infer that Ei and E2 are 
eigenvectors of Rp both with eigenvalue 1, since Pj = Rp. Appealing to the fact that Rp must 
satisfy the inequalities (21) and (22) above, since crp was assumed to be a scattering map, we 
deduce that Rifyp = -fp, whence the unit collision normal 'yp is another eigenvector of Rp with 
eigenvalue -1. 
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We restrict our attention to the subspace of 1.^ orthogonal to jp, namely 
Setting 

-(-^,0,-^,0,0,0) and £2 :-(o,-^,0,-^,0,o), (25) 

one may check that £1 • = E 2 ■'y /3 = 0, while evidently Ei ■ E 2 = 0. Let us consider any 

orthonormal basis of containing Ei and E 2 , namely {Ei,E 2 } U {ET,iP),EA,(fi),E^(fi)}, 

where each Ei{P) is allowed to depend on the spatial configuration /? e T^. One may then verify 
that any matrix of the form 

5 

Rp := ® + £2 ® £2 + ^ Ai(fi)Ei(fi) ® Eiifi) - 7/3 ® 7/3 

i=3 

with Ai(JS) € {-1,1) is a bijective linear involution which maps to Moreover, transforming 
back from quasi-momentum variables, the associated scattering matrix crp := M^^RpM conserves 
the total linear momentum and kinetic energy of its argument. The proof of the proposition follows. 

□ 


Evidently, as we have such a large family of scattering matrices which conserve both linear 
momentum and kinetic energy, it is prudent to specify another natural condition on each matrix crp 
to obtain a unique family of matrices {crp}p^j7, to which we can turn our attention. At this point, it 
is helpful to consider the case of hard disks. 


3.5.2. Comparison with the Case of Hard Disks. If we have developed a suitable extension of 
the classical scattering of hard disks to the more general compact, strictly-convex particle setting, 
the associated scattering matrix crp should reduce essentially to the classical Boltzmann scattering 
matrix (5) when P* is chosen to be a disk. We consider the case P* = B* (the closed unit disk 
in ]R^). As the classical Boltzmann scattering matrices are unique in the class of all maps on 
which conserve total linear momentum, angular momentum and kinetic energy of particles (and 
which also enforce non-penetration), we do not have mixing of pre-collisional linear velocities and 
angular speeds following collision. With this observation in mind, we consider the block scattering 
matrix defined on by 


' I4 - 2y{t//) ® 7 ((A) 

O2 ) 

, O4 

h ] 


g]^6x6 


for f e 


wifh 0„i, Im e R™^"* and 7(1/^) = -^[e(i/^), -e(i/^)]. Nofably, fhis mafrix is fhe identify map when re- 
sfricfed fo fhe sef E®. Mofivafed by fhis observation, we have fhe following corollary fo proposifion 
3.2 above. 


Corollary 3.3. Suppose crp is a scattering matrix satisfying the hypotheses of proposition 
3.2 which is the identity map when restricted to E^ = E^ n E^. Then crp is necessarily of the form 
crp = M~^{I - Typ ®'yp)M. 

Proof. Eel Bp = be any orfhonormal basis for E° which confains fhe vecfors E\ and 

E 2 given above by (25). Since by assumption crp\jp - l, if follows fhaf Ai(J3) = 1 for i = 3,4,5. 
Now, using fhe facl fhaf 
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we find Rp = I - ® 'yp. Transforming back to velocity variables V from quasi-momentum 

variables P, we obtain (Tp{V] = M~^ - 2yp ®'yp^ M, which yields the assertion of the corollary. 

□ 


As such, the derived family of scattering matrices reduces to the family of Boltzmann scat¬ 
tering matrices (which is the identity map when restricted to the factors of describing angular 
speed) when the reference particle P* is chosen to be a disk. With this concrete family of scattering 
matrices in hand, we now look to construct global-in-time classical solutions to Euler’s equations 
on D2. 


3.5.3. Construction of Global-in-time Classical Solutions on D 2 . We now offer some brief 
comments that establish Theorem 3.1, the proof of which follows swiftly from the construction of 
the scattering matrices crp - M~^{\ - 2yp i^'yp)M and an application of theorem 10 in Ballard 
[2]. We do not discuss technical details of the proof here, and refer the reader to ([2], section 
4) for details. Given that ()P* is of class C‘^ and that there is no externally-imposed force in the 
equations of motion (S^) and (S'*'), it follows that for each initial datum Zq e 2)2' there exists a 
unique piecewise linear map 1 1-^ {x{t),lc{t), P{t), Pf)] with 

[x(0),T(0),??(0),^(0)] =niZo and -^[x{t),x{t),mM)\ = n2Zo, 

at- t=o 

which satisfies (S~) and (S'*") on R and R \ 7” (Zq), respecfively. Moreover, for every such initial 
datum Zo € T'(Zo) the set of all collision times WlZo) is finite, i.e. TiZf) = with M = 

M{Zq) € N, with the property that for each t e (Ty,Ty+i], there exists a left-neighbourhood of 
t on which t i-> [x(t), x(t), ^^(t)] is analytic. Importantly, uniqueness of classical solutions 

allows us to define a hard particle flow {TdreR on 2)2 with the property that total linear momentum 
and kinetic energy of initial data is conserved for all time, and for which the colliding particles 
experience at most hnitely-many collisions on bounded time intervals. 

It is also important to emphasise that in order to make use of the general existence theory 
in [2], a family of scattering maps must be provided as data for the problem. As such, classical 
solutions are only unique with respect to the given family of scattering maps under consideration. It 
would be possible to construct another distinct hard particle flow on 2)2 that conserves total linear 
momentum and kinetic energy if one constructs a family of nonlinear scattering maps {crp}p^^7, on 
R® satisfying the same property. As intimated above, we do not address this problem in this article. 


3.5.4. An ‘Almost Physical’ Family of Matrices. It is important to record the fact here that 
the matrix up ■.= M~^ {l - lifp ® M € R®^^, where the unit vector€ R^ is given by 






-M" 


np 
-np 
Pp-np 
-q^ ■ np 


conserves the total linear momentum, angular momentum and kinetic energy of its argument, but 
it is not a scattering map in the sense of definition 3.7 above. In particular, one can find collision 
configurafions /?* e and associated pre-collisional velocities T* e 2^ \ 2° that satisfy 


up,V* = T* 


'To be precise, if Zq is taken to lie in dDi (namely the initial condition describes a collision configuration) then for 
consistency we should only allow for initial velocities n 2 Zo to lie in where/3 6 is determined by IIiZo. 
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and which therefore lead to interpenetration of the particles when the dynamics of (S“) and (S'*') is 
continued after collision. As a result, it cannot be used to construct a hard particle flow on 2)2 ^ but 
it can be used to construct a family of flow operators on corresponding to Euler’s equations 
presented in section 3.2. 

We make the rather naive comment that in the Boltzmann-Grad limit of the BBGKY hierar¬ 
chy, ‘particles become points’ and so it makes no sense to speak of non-penetration of particles for 
the limiting system as the number of particles A —> 00. As such, one could argue that the family of 
maps would nevertheless be suitable to establish a kinetic model for the average behaviour 

of rarified gases composed of compact, strictly-convex particles. Indeed, the main result Theorem 
4.1 on characterisation of collision invariants for non-spherical particles in this article also holds 
for the family under the weaker condition that 5P* be of class C\ as opposed to analytic. 


4. Collision Invariants for Compact, Strictly-convex Particles 


We now turn to the proof of the main result of this article. We firstly define the analogue of 
classical collision invariants in the case when the underlying particles are not disks. 

Definition 4.1. Let S = be a family of maps on R®. A measurable function (f : 

X R X — > R is said to be an S-collision invariant if and only if it satisfies the functional 

equation 

a)'p, 1?) -I- aj'p, ’d) = (p{v, oj, §) + tp(y, To, (26) 

for every V = [v, v, co, ui] e R^, = {§, ijj) where 


crpiV], \ ^ / crpWh \ 


We also make one more definition. 


Up := crpiV]^, Up := crp[V](,. 


Definition 4.2. We define to be the class of reference particles P* c R^ which have 

reflection symmetries in the two canonical orthogonal axes of R^. 


We are now ready to state in precise terms the main result of this article. 

Theorem 4.1 (Characterisation of Collision Invariants). Suppose P* € has the property 

that 5P* is analytic, i.e. of class C^. Let S be the associated family of matrices 

-2yp®yp)M}p^^3 or -2 pp®'7 ]p)M}p^^},. 

If a measurable map is an S-collision invariant, then it is necessarily of the form 

(p{v, u, d) - a{b) -\-b ■ V + c (m|vp -I- Ju^^ , 

for some , Z12, c e R and some measurable a : —> R. 


Remark 4.1. For the proof of this theorem, we need only restrict our attention to the family of 
maps {cTgl^gx^ with crp = M~^{l-2flp®yp)M, as the proof for the family crp = M~^{I-2j]p®'rjp)M 
follows directly thereafter. This will become clear in section 4.4 below. 

4.1. Rewriting the Functional Identity. It will prove useful to rewrite the functional iden¬ 
tity (26) in a way that allows us to employ properties of reflection matrices, to which each scattering 
matrix crp - M~^{1 - 2y^ ®'yp)M is conjugate. Indeed, given a collision invariant ip, we define 


^o(v, u,-&) (p{v, u, d) - (fiO, 0, d), 
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together with an associated map x —> Rby 

V; ?9-, ?9') := tpo{v, co, &) + (po(y, To, &), 

noting that d>;p(0; ■&,■&)- 0 for all {■&, &) e T^. It follows that ^ is a collision invariant if and only if 

%{(Ti 3V\ I) = %iV-, I) (27) 

for all yS e T^. Setting P ■.= MV, and also define the new map O* : R^ x —> R by 

I) ■- 

we find fhaf i/j is a collision invarianf if and only if 

(28) 

for all /? € and P e R®. If is now we make fhe imporfanf observafion fhaf if fhe orienfafions 
{•&, •&) we fixed, fhen (28) implies fhaf 

d)*(-; P, ■&) is consfanf on fhe group orbifs G^P, 

for any chosen P € R^, where c 0(6) is fhe subgroup generafed by fhe 1-paramefer family of 
reflection mafrices {/ - ly^ : \fj e S^), namely 

Gf-({/-2 y^®r/y : 

Transforming back fo T-variables and observing idenfify (27), we obfain fhe following resulf. 

Proposition 4.2. The map ip is a collision invariant if and only if for each (•&, &) € T^, the map 
<1>^(-; P, P) is constant on the group orbits Q^V c Bfi for every V e R^, where 

The basic problem is now fo characferise fhe orbifs of every poinf in R^ under fhe acfion of 
G^ for each (P, P) € T^. This leads us fo fhe concepf of energy-momentum submanifolds of R^, 
which we infroduce now. 


4.2. Energy-momentum Submanifolds of R^. We recall from secfion 3.4 above fhaf fhe 
scaffering mafrices o-p under sfudy conserve fofal kinefic energy, i.e. 

IMo-pVi^ - |MPp, 


along wifh fofal linear momenfum of fhe particles, 


m 


(^fiV)i 

{CTpV)2 


m 


((rpV)3 

(CTpV), 


= m 


El 

E2 


-I- m 


y^ 

Va 


for all yS € T^, once P € R^ has been prescribed. Writing fhese in fhe language of fhe previous 
secfion, we have fhaf 


Y 1-^ |MTp is consfanf on fhe group orbifs Q^V 


fogefher wifh 

T i-> Ti -I- F3 and T i-> T2 + ^4 are consfanf on fhe group orbifs Q^V, 
for T € R®. On fhe basis of fhese observafions for fhe family of scaffering mafrices 

If/p ® 7y3)M)^gT3- 

if seems reasonable fo posfulafe fhaf fhe orbifs Q'^V are simply fhose subsefs of R^ which are 
realised as fhe infersecfion of energy ellipsoids 

&(V) ■- (f € R^ : |MFp - \MV\^] 
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with momentum planes 

ni(l/) jy e : yi + Fa = Vi + 1/3) and HaCF) {f € R^ : F2 + F4 = F2 + F4) . 

Indeed, this is what we prove in section 4.4 below by using group-theoretic techniques and a careful 
analysis of properties of the collision normals jp. 

Firstly, let us put the previous argument in precise terms. For a given energy e > 0 and 
momentum vector p e R^ satisfying > |pp/2m, we define the associated energy-momentum 
manifold M(e, p) c R^ by 


M(e, p) F € R® : \MY\ = e and 


i^i + Fa 

F2 + F4 

It is now our aim to show that if F € R'’ is arbitrary, and we denote 


then the group orbits of points F € R^ are given by 


giv = 


M(e,p) 

l[tt0.o]) if8^ = S 




P 

m 


for any choice of orientations {■&, d) € T^, i.e. the group orbits are independent of the choice 
of particle orientations. In other words, we want to show the restriction map d)(-; i9')lM(e,p) is a 

constant function for all suitable e > 0 and p e R^ by identity (27). Since we have nothing to 
show in the case that is a singleton set, we assume henceforth that > |pp/2m. As the 
energy-momentum submanifolds are homeomorphic to the 3-sphere, one can expect to reduce the 
study of the subgroup c 0(6) acting on M(e, p) to one of a group acting on S^. As done before 
in section 2.3 above, let us now reduce our problem to a kind of canonical form. 


4.3. Transformation to Canonical Form. Let both energy e and momentum p be given 
which satisfy > |pp/2m, and suppose them to he fixed. We now define /ie,p : M(e, p) —> by 




1 

W) 


(MF)i - (MF )3 ' 

(MF)2 - (MF)4 

V2(MF)5 

V2(MF)6 . 


where r(F) := ^((MF)i - + ((MF)2 - {MV) 4 f + 2{MV)l + 2(MF)2, thereby consider¬ 

ing as embedded in R^. Notice also that since > |pp/2m, the radicand of r(F) is strictly 
positive. One can check that h^ p is a bijection between M(e, p) and S^, whose inverse is given 
explicitly by 


k!pM = j 


J2e2 - ^wi + ^ 

\ fn ^ ^/m 

4 


yfm 

W2 + 

ni ^ ^Jm 


2e2-^ 

m 

Pi. _ ./2e2 - 
V7i ^ 

^ _ /2e2 - ^ 

V 2 4m 2 


Wi 


>V2 


for W = (wi, W2, W3, W4) € 
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We now consider the orbits Q^V as images of another group action on under the map hg p. A 
quick calculation reveals that 

(Tp ■- (/ - Ijp ®yp) MV for V € M(e, p) 


if and only if 


sp : w i-> ^7 - l/up ® fip^ w for w = he^p(V), 


where np € is the unit vector 



{Np)i 


(Np)2 

^ [rp - dpe{iP)f ■ Np 



A/J > 0 is given in (18) above, and the reduced mass-inertia matrix M\ € is given by 


Ml 


^Jm 0 0 O' 

0 ^fm 0 0 

0 0 V7 0 

0 0 0 V7, 


It will be crucial for the proof of characterisation of collision invariants in the sequel to show 
that the (?7, j^l-dependent family of unit vectors (pp : € S^) lies in no single hyperplane in 

R^. Indeed, we address this problem in proposition 4.5 below. With this observation that we may 
essentially work on the sphere for any pair of orientations {&, d) e T^, we define the group 
hI c 0(4) by 

hI ■.= {[!-iPp^'Pp : (A€S')), 

which is now the primary object of study. We have the following proposition, which crystalises 
the above discussion. 

Proposition 4.3. Let (i7, ■&) e be given. The group 77^ c 0(4) acts transitively on if and 

only ifQ^ c 0(6) acts transitively on M(e, p)/or any single pair (e, p) satisfying > |pp/2m. 

If the orbits under of any given point in R^ is indeed the corresponding energy-momentum 
manifold, we may immediately infer the existence of another measurable function : R^xR —> R 
such that 

^(p{V;&, d) = (^ip{mv + mv, m\v\^ + Joj^ + mfv'^ + Jaj^), (29) 

for all P e R®. To show that (29) holds for some <1)^, we employ some new results contained in the 
appendix of this article on generators of the rotation group 0(4), which are due to C. Viterbo. 

4.4. The Transitive Group Action of 77^ on The key result is the following, whose 
proof can be found in Appendix A. 

Theorem 4.4. Suppose that p : is a continuous, non-constant map. Let 77 c 0(4) 

denote the group 

77 ^|7 - 2p^ f . 

Then 77 acts transitively on unless the image set [p^ : e SM A strictly contained in some 

hyperplane in R^. 
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Using this result directly, we are able to prove that defined above does indeed act transi¬ 
tively on In fact, the proof of Proposition 4.3 follows immediately from the following result, 
which says that the image set (pp : i/r € S') cannot lie in any one fixed hyperplane for any choice 
of orienfafions {■&, i?) € T^. 

Proposition 4.5. For any {&, F) e T^, we have spanl/T^j : i/^ e S') = 

Proof. Lef {■&, &) e be given. We suppose, for a confradicfion, fhat fhere exisfs a vector 
W{&, ■&) \ (0), written componentwise as 

W{&, - (wi(i9', F), W2{'&, &), W3(i9', F), wpF, F)), 


such that 

Pp-W{F,^) = 0 for all (A € s'. 
However, we note that this is equivalent to the statement that 

yp-V{-&,d-) = Q for all (A e s'. 


(30) 


where 

¥{■&,■&) := yJTJmwi,- yjlimw2, pljJwj,, . 

Importantly, assumption (30) implies that 

n 


where n In particular, there is at least one non-zero velocity vector Vph, &) which is 

both pre- and post-collisional for every choice of elevation angle (A e S'. As the notion of pre- 
and post-collisional velocities is inherently dynamic, we must now appeal to the existence results 
established in section 3.5. 

Let us consider the following 1-parameter family of initial data Zo(iA) = bo>Zo(<A)] (parame- 
terised by lA € S') for the ODEs presented in section 3.2, where 

Zo = [0, ■&, 0, oj] and zo(iA) = [df{il/)eiil/), v, oi]. 


with 

“•* “^=-V7"’4- 

It then follows that for the initial data [zo>Zo(iA)] their associated phase trajectories t i-> x^{t), 
t 1 -^ and t 1 -^ Icpt), t i-> i9'i/,(f) (which are smooth, by the results in [2]) there exists d > 0 
independent of (A such that 


F{x^{t), Xi/,(t), Fifrit)) > 0 for all - d < t < d. (31) 


We use this deduction to derive our contradiction by reducing our considerations to properties of 
the motion of the point of contact on particle P both before and after collision. In the sequel, we 
often suppress the dependence on &, ■& for all relevant quantities of interest, in order to make the 
presentation of our arguments clearer. 

We perform a time-dependent change of variables so that particle P is stationary for all time, 
and the dynamics of P takes place in the exterior domain \ P. It will be convenient to take 
the view of material point trajectories which evolve on the particles P and P. Firstly, let Xp{f, xf) 
denote the position of the point on particle P at time t e R whose initial position at time t = 0 is 
xo, namely 


Xp{f,xo) R{ajt)xo 
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for any xq e P- Similarly, let Xp(t; xq) denote the position of the analogous point on particle P at 
time t € R, i.e. 

Xp ( t ; xq) := R{ojt) (xq - d{i]j)e{i]j)) + d{il/)e{il/) + vt, 

for any xq e P. Transforming to the time-dependent reference frame from which P is viewed as 
stationary, Xp(f; xq) Xp{t, xq) and Xp(t; xq) X-^t; xq), where 

Xp{t,Xo) - Xq for Xq G /?(i9')P*, 


for all time t e R, and 


Xp(t; a:o) = R{{co - co)t) (xq - d{itj)e{\jj)) + R{-ajt) {d{ip)e{ip) + vt ), 

for Xq € /?(?9'(t))P* -i-v(t). As such, we may conveniently view the motion of individual points on the 
particle P as taking place in the exterior domain R^ \ /?(??) P*. In order to derive our contradiction, 
namely that must indeed be the singleton {0}, we focus our attention on the trajectory of 

the point of collision which lies on particle P(t). For the C^(-d, d) trajectory t i-> Xp{f, p{t)/)) to 
satisfy 

{%(t; p((A)) : t e (-d, d)) c R2 \ /?(j?)P* for all lA € 

it is necessary that the normal component of the curve {Xp{t', p{^)) : t e (-d, d)) vanish at f = 0, 
i.e. 

—Xp(t; p{ifj)) ■ n{ij/) = 0 for all lA € 
dt ^ t=o 

A calculation reveals that this holds if and only if 


^ii//)-W = 0 foralllAeS^ 

where ^ = ^^(lA) £ is given by 



_ «^(<A) _ 

(p^(iA) - ^/^(lAMiA)) -w^CtA) _ 


(32) 


We now show that the linear span of the set |^^(iA) : lA ^ is the whole space R^ for any choice 
of (d, d) € T^, which implies that IT € R^ must indeed be the zero vector by (32) above. We 
require the result of the following simple lemma. 


Lemma 4.6. Suppose P* e There exist at least two angles lAi - i]j\{d,d),i]j 2 - 

^ 2 {d,d) e such that p^iihi)-^ ■ n^{i//i) = 0. 


Proof. We recall that one axis of symmetry of P* lies along the x-axis, and the other lies 
along the y-axis. We denote by djc > 0 and d^ > 0 the largest positive values of the x- and y-co- 
ordinates that lie on these axes of symmetry, respectively. Consider the angle lAi = if/i{d,d) e 
that gives rise to the point p^(iAi) - R{d){6x,0) and the associated normal vector n^(iAi) to P at 
p^ftpi). Since the reference particle P* has Z 2 x Z 2 symmetry, it follows that/?(?9').S'i/?(i^)^P = P. 
Moreover, as dP* is of class C" and so the outward normal at p^(iAi) is unique, it follows that 
^^(<Ai) = ^(d)(l,0), whence p^CtAi)"*" • n^(iAi) = 0- Ths other case follows by considering (A 2 = 
d + njl, and arguing similarly by using the fact that R{d)K 2 R{d)^P = P. □ 


We now make the following four judicious choices of the angle of elevation lA e to produce 
vectors {^ 1 ,^ 25 ^ 35 ^ 4 } which are candidates for a basis. Using the result of the above lemma, we 
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choose e with the property that - R{&){6x,0) and = 0, which yields 


the vector ^?(iAi) given by 


^i = Qi? 


Qi? := 


where € 0(4) is the rotation matrix 


1 

0 

0 

r/^((Ai)sintAi 


^ cos - sin I? 0 0 
sin I? cos I? 0 0 


0 

0 


0 

0 


1 0 
0 1 


Choosing 1/^2 = *Ai + ^/2 and following similar reasoning, we yield ^2 - ^{^ 2 ) given by 


^2 - Qi5 


0 

1 

0 

dlifPi + f)siniAi 


Next, we choose any with tpi < tjjj, < 1^2 satisfying the property that 


_r/^((Ai)sin^i 
dfMi + I) sin ij/i 




together with ^ 0 ; we note that this is always possible since p^{i//)'‘~ ■ n^(iA) = 0 

for all ij/ satisfying < tp < ^2 + ^/2 if and only if P* is a disk. Indeed, for such a 1/^3 e S^, we 


set ^3 Q{p 3 ), where 


^3 Qi? 


- d)2 

■ «^(<A3) _ 

(;^^(<A3) - d^{ip3)e{ip3)) ■ n^{il/3) 


Finally, we choose >p^ = p 3 +n and set ^4 = APa), which yields by symmetry that 


^4 Q & 


-«^_^(<A3 - ^?)l 
-«^_^(</'3 -d)2 

■ «^(<A3) _ 

(f^(<A3) - d^{ip3)e{il/3)) ■ n^{il/3) 

With these observations in place, we approach the following lemma. 

Lemma 4.7. The set o. basis for if and only if sin ipi 0. 

Proof We need only show that (4^',^',^',^ 4 ) is a basis for R^ when sini/^i 3 ^ 0, where := 
Qj^y. Evidently, is a set of linearly independent vectors. Assume for the moment there 

exist constants (ci,C 2 ,C 3 ) € R^ \ { 0 } such that 


^4 - ^ 1^1 + ^24^2 + <^3^3- 
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By necessity, C 3 = 1, since 1/^3 € S' was chosen so that + 0. This immediately 

yields that ci = - 2 n^_^(i /^3 - ■&)i and C 2 = -2n-^_^{i]j3 - '&) 2 - However, with these values of 
constants c, it must be that 


_d^(iAi)sini/ri 

d^(^\ + f)sin(Ai 


-d) = 0 , 


which contradicts the properties of the elevation angle 1 /^ 3 . Thus, cannot be a linear combination 
of , ^' 2 , and so the set {^'constitutes a basis for in the case where sin tj/i 0. □ 


To conclude the proof of the proposition, we need to consider the construction of another 
basis in the case when sini/^i = 0. To do this, we consider the auxiliary function on §' given by 
the rule 

If/ I—^ . 

Notably, this function vanishes when i// = i/zi or i// = if/ 2 . Importantly, the numerator and denomi¬ 
nator are both identically zero for all i/r € S' if and only if P* is a disk. Since, by assumption, P* is 
not a disk and its boundary 5P* is C", this function is non-constant and smooth away from those 
points where the denominator vanishes. We therefore choose any two distinct 1 /^ 3 , 1/^4 € S' with the 
property that p^iif/i)^ ■ n^{tj/i) 0 for / = 3,4 and 

• nj_^iif/3) d:^_^iif/4)e(tf/A)-^ • n:^_^iif/4)' 

Using this observation, it follows from an argument identical to that found in the proof of lemma 
4.7 that the family {^\,^ 2 ,^ 3 ,^] constitutes a basis for 1.^, where 





«^_,?(‘A4)i 

^3 ■= Qi? 


0 

Jl 



Thus, we have shown that the span of the set {ppii]/) : i/^ € S'} is indeed R^, which completes the 
proof of the proposition. □ 

We conclude by noticing that by Theorem 4.4 the group acts transitively on S^ for every 
6 e s', which immediately yields that acts transitively on energy momentum manifolds for any 
choice of orientation pair j?) e T^. As a result, there exists a measurable map such that 

■&) = <l>y(mv -I- mv, m|v|^ -I- + m|v|^ -I- Tio^). 

We now prove that this representation formula implies that collision invariants ip are necessarily of 
the form 

p{y, a>, d) = a{&) + b -v + c ^m|v|^ + , 

for any constants bi,b 2 ,c € R and any function a : S' —> R. To do this, we appeal to classical 
results on Cauchy’s functional equation. 


Remark 4.2. We believe that proposition 4.5 holds true for an arbitrary compact, strictly 
convex reference particle P* in R^ with C" boundary, although we have chosen not to explore this 
particular extension of proposition 4.5. 


Remark 4.3. As one need not appeal to dynamical considerations in this case, the proof of 
proposition 4.5 also holds for the family of matrices {M“'(/ - Irfp when the boundary 

curve 5P* of the associated reference particle 5P* is only of class C', as opposed to analytic. 
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4.5. Cauchy’s Functional Equation. The last remaining step in the proof of Theorem 4.1 
is proving the following proposition. 

Proposition 4.8. Let e > 0 and p € such that > |pp/2m, and let ip be a collision 

invariant. Suppose that i9')lM(e,p) ^ constant function. Then p is necessarily of the form 

p{v, CO, d) = aid) + b -v + c {in\v'^ + Ja/^ for V € R^, 
for constants Zii, Z 12 , c € R and a function of orientation a : S' —> R. 

Proof. The main idea of the proof is to transform identity (26) for collision invariants into 
Cauchy’s well-known functional equation for a real-valued function g on R, namely 

g(x) + g(y) = g(x + y) for V, y € R. (33) 

It is well known (see Darboux [7]) that under the assumption g he continuous at a single point of 
R, the only possible solutions of (33) are linear functions g(z) = cz, where c e R. Since we assume 
p to be measurable, Lusin’s theorem immediately gives us enough continuity of (/> on R^ x S' for 
the following arguments to be valid. Indeed, since (26) is equivalent to identity (27), we notice 
that if d>i^(-; d, ?9^)lM(e,p) is constant then is necessarily of the form 

^ip{V\d, d) - 0^(mv + mv, m|vp + Jo/ + m\vf + laf", d, d) 
for some measurable auxiliary function d)^. Since it then holds by definition of d)^ that 
d)y(mv -I- mv, m\vf + Jo?' + m|vp + Jo?', d, d) - po{v, 00 , d) + po(y, m, d), 
setting V = 0 and ai = 0, we find that 


poiv, 0 ), d) = d)^(mv, m|vf + Jo)^\ d, d), 

namely that the value of d)^ is independent of its second parameter d. By repeating this argument 
by instead setting v = 0 and a» = 0, we conclude that d)^ is independent of both d and d, namely 
that 

Pq{v, 0), d) + pQ(fi, Td, d) = 'Tp(v -I- v, \vf + ^ 0 ? + |vp -I- jfo?) (34) 

for some new measurable function Since p is assumed to be a collision invariant, if follows 
that satisfies the identity 

|vp -P ^m^) -p ^^(T, |vp -p jfU?) - 'P^(v -P V, |vp -P j^o? -P |vp -P jfd?). 

Finally, setting a> = ai = 0, we infer that 

'T^(v, |vp) + 'T^(T, |Fp) = ^^(v + iJ, |vp + |vp). (35) 

It is at this point we invoke an argument from Truesdell and Muncaster [15]. Let us now make 
the choice v = -v, which yields from (35) that 

^-^(0,2|vp) = ^^(v, |vp) + ^^(-v, |vp). (36) 

Next, selecting any two orthogonal vectors v, v, we deduce from (35) that 

T„(0,2|vp + 2|Tp) = ^.(0,2|v + vp) 


(36) 


'T„(v + V, |vp + FP) + T„(-v - T 


TT l,,|2 , it:i2 


+ 




'P^(v, |vp) + T'^(v, |vp) + T'^(-v, |vp) + 'T^(-v, |vp) 
'P^(0,2|vp) + T'^(0,2|vp). 


( 37 ) 
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Thus, the map gi(i') T'^p(0, s) satisfies Cauchy’s functional equation on [0, oo), and is therefore 
necessarily of the form gi(i') = cs for some c e R. Now consider the map g 2 (v) T^(v, |vp) - 

gi(|vp). One may check that g 2 is measurable and odd on and by (35) and (37) above is additive 
on orthogonal pairs of vectors in R^. It follows from (Truesdell and Muncaster [15], page 88) 
that g 2 is necessarily of the form g 2 (v) - b ■ v for some b e R^. As T'^(v, |vp) = gi(|vp) + g 2 (v), it 
follows that 

|v|^) = b-v + c|vp. 

Thus, setting v - 0 and ai = 0 in (34) above, we deduce that (fio satisfies 

99o(v, CO, &) - b ■ V + c + Jco^^ 

for some b e and c € R. Since any funcfion of ?? € S' is a collision invarianf, fhe claim of fhe 
proposifion is proved. □ 
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Appendix A. On Groups Generated by Reflections (by Claude Viterbo) 

We shall here prove fhe following resulf abouf fhe fransifive group acfion of H'J on S^. Lef 
yu : s' —> S^ be a continuous curve and 5 : S' —> 0(4) be fhe associafed hyperplane symmefries 
wifh respecf fo jU"*", namely - I - l/i^p ® /i,/, for i/r e S'. 

Proposition A.l. The group generated by the reflection matrices {s^ : i/i e S') acts transi- 
tiveiy on S^ uniess the image of p is contained in a hyperpiane ofRfl. 

If is imporfanf fo menfion fhaf proposifion A.l exfends fhe work of Eaton and Perlman ([8], 
theorem 1), in the sense that we do not need to take the Euclidean closure of ({5,^ : e S')) 

in order to infer that it is indeed the whole group 0(4). In what follows, we actually prove the 
following more general result, from which A.l follows. 

Proposition A.2. Let p : A ^ S”“' be a continuous map, where A is connected and not 
reduced to a point. Let s^ := I - 2p^ ® p^|, be hyperpiane symmetry matrices with respect to p-^. 
The group G generated by {s^ : i/' € A) c 0{n) is identicaiiy equai to 0{n) uniess there is a 
k-dimensionai hyperpiane ncR”(k<n-l) such that p^ € If for aii if € A 

Note that if the image of p is contained in a hyperpiane H, the orthogonal set H-*- is invariant 
by all elements of the group generated by [ 5 ,/, : i/i e A) c 0{n) and thus the associated action on 
S”“' cannot be transitive. 

The following result has been proved in [8] : if the group G is infinite, then its closure is equal 
to 0(n). But since the map s is non constant, the group generated by the elements s^ is necessarily 
infinite. We may thus assume G is dense in 0(n). 

We note that the hyperpiane symmetries s^ have determinant -1. It will be useful to consider 
the group K, the intersection of G with all proper rotations of 4-space SO(4). Since every element 
of K can be written as the product of an even number of matrices s^, we have the following: 

Eemma a.3. The group K is arcwise connected 
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Proof. As SMs arcwise connected, we have g = homotopic to = 

= / in SO(4) for any g & K, where 1 denotes the identity element of □ 

We shall also need the following theorem. 

Theorem A.4 (Kuranishi-Yamabe-Goto). Let H be any connected subgroup of a Lie group 
G. Then H is a Lie group. Moreover, there is a Lie subalgebra h o/ 9 such that there exists a 
neighbourhood V of the identity e in H and f/ c h with V - exp(h n U) 

Proof. We refer the reader to [10] (see also [9] Theorem lip. 292, and p. 196). □ 

Finally if G is a Lie group and H a connected subgroup, there is a (proper) maximal connected 
subgroup of G containing H. We do not require Zorn’s lemma, since we may simply take a 
subgroup of maximal dimension strictly less than dim(G) containing H. 

Lemma A.5. A connected maximal subgroup ofSO{n) is necessarily closed, hence compact. 

Proof. A maximal subgroup is either closed or dense. We could use ([14] Theorem 1.3, p. 
628) applied to the special case of SO(n), which has the property the the connected component of 
its center is trivial. This result states the following: if G is a connected Lie group and /i : G —> 
SO(n) is a Lie group homomorphism with dense image, then h{G) = SO{n). □ 

Remark A.l. In the case n 4 when the group SO{n) is simple, we have a simpler proof. 
Indeed, according to Theorem A.4, such a subgroup corresponds to a Lie algebra of so(n). Let 
then h be a Lie subalgebra of so(n) corresponding to a dense subgroup H. Since Ad(g)I) = h for 
all g e //, we have by density that this still holds for any g € SO(n), hence h is an ideal of so(n) 
and // is a connected normal subgroup of SO(?i). But this is impossible, since so(n) is a simple Lie 
algebra. 

Proof of Proposition A.I. The group K is dense, connected, and contained in a maximal con¬ 
nected subgroup which is of course dense. Thus K = SO{n). It is then follows at once that 
G = 0(n). □ 
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